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ABSTRACT 


The  velocity  distribution  of  a  point  excited, 
infinitely  long,  thin  cylindrical  shell  immersed  in 
an  acoustic  medium  is  considered.  The  problem  is 
analyzed  by  applying  integral  transform  techniques; 
a  solution  is  then  obtained  by  evaluating  the  Inverse 
transform  Integrals  numerically.  Examples  are  presented 
to  Illustrate  the  effect  of  fluid  loading  on  the  vibra¬ 
tional  behavior  of  a  point  excited  cylindrical  shell. 

In  addition,  results  are  compared  with  those  of  a  point 
excited  plate  with  and  without  fluid  loading. 


ADMINISTRATIVE  INFORMATION 

The  Investigation  presented  in  this  report  was  initiated  under  the 
Independent  Research  Program,  Program  Element  61152N,  Task  Area  ZR0110801, 
Work  Unit  1960-011  and  completed  under  funding  from  the  Naval  Sea  Systems 
Command  (OSH)  under  Task  Area  SF4 3452702,  Task  18185,  Work  Units  1960-010 
and  1960-020.  ~ . . .  — 


INTRODUCTION 

Understanding  the  vibration  and  radiation  characteristics  of  a 
cylindrical  shell  is  of  considerable  importance  to  the  Navy.  It  can  provide 
insight  Into  the  reasonableness  of  simplifying  approximations  currently 
utilized,  for  Instance,  determining  the  frequency  range  in  which  a  complex 
structure  can  be  modeled  as  a  flat  plate.  Unfortunately,  the  analysis  of 
a  fluid  loaded  cylinder  is  complicated  by  the  fact  that  there  are  three 
components  of  displacement  which  are  coupled  to  each  other.  This  differs 
from  a  flat  plate  where  only  one  component  of  structural  displacement  is 
considered. 
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Determining  a  solution  to  the  shell  equations  for  an  Infinite  cylinder 
is  further  complicated  when  fluid  loading  is  included.  However,  the 
response  of  an  infinite  cylindrical  shell  with  or  without  fluid  loading  can 
be  determined,  in  theory,  by  applying  Integral  transform  techniques.  The 
difficulty  with  this  method  is  that  complicated  Inverse  transforms  must  be 
evaluated.  Past  studies  have  circumvented  this  difficulty  by  several  ap¬ 
proaches.  One  is  to  consider  the  exciting  force  specification  to  be  inde¬ 
pendent  of  the  coordinate  defining  the  axial  direction  of  the  cylinder, 

1* 

which  leads  to  a  plane  or  two  dimensional  problem.  Other  investigators 

have  simplified  the  problem  by  assuming  an  exciting  force  which  is  a  periodic 

2  3 

function  having  a  specified  axial  wave  length.  ’  However,  once  a  longi¬ 
tudinal  wave  length  is  specified,  then  only  one  longitudinal  wave  length 
is  permitted  for  the  response  of  the  infinite  cylinder.  But  in  reality, 
all  wave  lengths  are  permissible  for  a  point  driven  infinitely  long 
cylindrical  shell.  Another  configuration  which  avoids  solving  the  inverse 

transform  integrals  is  that  of  a  cylindrical  shell  with  periodic  stiffeners 

4 

which  reduces  the  Integral  to  a  summation.  In  any  event,  utilizing 
any  of  the  above  simplifying  procedures  negates  to  a  varying  degree 
the  usefulness  of  the  solution  to  the  actual  physical  problem.  On  the 
other  hand,  with  the  use  of  large  computers  and  refined  numerical  techniques, 
it  is  now  possible  to  determine  the  response  of  a  point  excited  fluid  loaded 
cylinder  or  its  radiated  pressure  directly,  thus  avoiding  the  above  inhib¬ 
iting  assumptions. 


*A  complete  listing  of  references  is  given  on  page  53. 


FORMULATION  OF  PROBLEM 


DESCRIPTION  OF  PROBLEM 


The  problem  considered  here  consists  of  an  infinitely  long  cylindrical 
shell  with  a  radius  "a"  much  larger  than  the  shell  wall  thickness  h.  The 
cylinder  is  surrounded  by  an  infinite  fluid  medium  with  a  fluid  density  p 
and  sound  speed  c.  The  cylinder  is  excited  by  a  temporally  varying  point 
force  Fe  where  F  is  the  amplitude  of  the  force  and  U)  is  the  circular 
frequency  of  excitation.  It  will  be  assumed  that  the  cylinder  is  in  steady 
state  vibration  and  therefore  a  time  dependence  of  the  form  e  for  all 
dependent  variables  is  assumed  throughout  the  analysis. 

The  cylinder  is  defined  by  an  r,  ()>,  and  z  coordinate  system  where  the 
r  axis  is  directed  outward  along  the  normal  to  the  shell's  surface,  (}) 
defines  the  angle  in  the  circumferential  direction  and  z  defines  the  longi¬ 
tudinal  axis  as  shown  in  Figure  1.  The  above  orthogonal  displacements  of 
the  shell's  middle  surface  are  denoted  by  w,  v  and  u,  the  radial,  tangential, 
and  longitudinal  (or  axial),  respectively. 


SHELL  EQUATIONS  OF  MOTION 

For  a  thin  cylindrical  shell  the  equations  of  motion  are 


3  u  ,  1-v  3^u  ,  1+v  3^v  .  V  3w  ii 

2  2  2  ^  +  -„-U 
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where 


,  the  compresslonal  wave  velocity  in  an 
elastic  plate 


P^d-v-) 


E  =  Young's  modulus 
V  =  Poisson's  ratio 
D  =  mass  density  of  shell  material 


h  =  shell  wall  thickness 
p^  =  total  pressure  acting  on  shell 
a  =  radius  of  cylindrical  shell 
The  total  pressure  acting  on  the  shell  can  be  expressed  as 


(Ic) 


F5(z)6(([))e 


■lu)t 


-  p(a,(|>,z,t) 


(2) 


where 

p(a,((),z,t)  =  p  (a,(ti,z)e“^“*^ 
o 

6((j))  =  delta  function  with  argument  4) 

The  first  term  on  right  hand  side  of  Equation  (2)  is  due  to  the 
harmonically  varying  point  force  whereas  the  second  term  is  due  to  the 
fluid  pressure  acting  on  the  shell  surface. 
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The  motion  of  the  shell  can  be  described  by  assuming  displacements 
of  the  form 


u(4),z,t) 


oo 

n=0 


cos  nd)  U  (z)e 
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-lojt 


v(4>,Z,t) 


w((t>,z,t) 


CD 

=2:- 


sin  n(p  V  (z)e 
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-iwt 


n=0 


uu 

=2’< 

n=0 


cos  nd)  W  (z)e 
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-iLi)t 


(3a) 


(3b) 


(3c) 


FOURIER  TRANSFORM  METHOD  OF  SOLUTION 

Solutions  to  problems  of  the  type  considered  here,  that  Is,  systems 
described  by  independent  variables  extending  over  an  infinite  domain, 
are  greatly  facilitated  by  employing  integral  transform  techniques. 

The  Fourier  transform  pair  which  will  be  used  in  this  analysis  is  defined 


F(k) 


+00 


'^00 


(4) 


which  transforms  the  f(z)  into  the  wave  number  domain,  and  the  inverse 
transform,  ^ 


f  (z) 


^  J  e^’"‘'F(k)dk 


(5) 


which  transforms  the  function  F(k)  back  into  the  spatial  domain. 

Substituting  Equation  (3)  into  Equation  (1)  and  transforming  the 
result  into  the  wave  number  domain  yields 
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U  (k)  -  ikan 


(1+v) 
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V  (k)  -  ikavW  (k) 

n  n 
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(6a) 
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(1+v)  .  r,  2  2  (1-v)  .  2  ^2 

ikan  ^  U  (k)  +  k  a  ^  +  n  - 

^  n  L  2 


V  (k)  +  nW  (k)  =  0  (6b) 

n  n 


IkavU  (k)  +  nV  (k)  + 
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2ir 


9  2  2  2  2 
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cos  n4»d(j)  (6c) 


where  ^  is  a  nondimensionalized  frequency  parameter. 
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The  above  equations  were  obtained  by  utilizing  the  orthogonality 
conditions 

2it 
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mn 


2TT 
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(7a) 
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where 


1  n=m 

6  =  {  n  I  and  E  = 

nm  1  0  nfm  n 


1  n=0 

2  n>l 


The  term  F6((}))  on  the  right  hand  side  of  Equation  (6c)  due  to  the 
exciting  force  can  be  expanded  in  terms  of  a  Fourier  series,  i.e. 


F6(({))  cos  n(J) 

n*0 


(8) 
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Multiplying  both  sides  by  cos  mij>d4>  and  integrating  between  0  and  2Tf  yields 

2ir 

F  =  cos  m(J)d<|)  “  F  (9) 

0 

Consequently,  from  Equations  (8)  and  (9) 

00 

F6((|))  cos  n(J)  (10) 

n=0 

PRESSURE  DUE  TO  FLUID  MEDIUM 

The  term  P(a,(}>,k)  in  Equation  (6c)  is  the  Fourier  transform  of  the 
pressure  acting  on  the  shell  due  to  the  presence  of  the  fluid  medium.  This 
term  can  be  related  to  the  radial  motion  of  the  shell  by  utilizing  the 
wave  equation  governing  the  fluid  medium  and  the  appropriate  boundary 
condition  at  the  fluid-shell  interface.  The  wave  equation  in  cylindrical 
coordinates  is 


,2  .  ,  , ,,  19  p{r.(J),z.t)  - 

1  p(r,<t),z,t) - ~  =  0 

c  at'^ 


(11) 


where 
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Suppressing  the  harmonic  time  dependence  the  wave  equation  can  be 


written  as 


V^p  (r,(|),2)  +  k^p  (r,()),z)  =  0 
o  o  o 


where  k  =  —  Is  the  wave  number  of  the  acoustic  medium, 
o  c 

To  determine  the  transformed  pressure  P  (a,<J);k)  the  z  domain  In  the 

o 

wave  equation  Is  transformed  Into  the  wave  number  domain. 


/  3^  13  13^  1,2  ,,2\  /-  A 

V3r^  r  3r  r^  30^°  /  ° 


k)  =  0 


Since  the  fluid  medium  Is  Infinite  in  extent,  there  are  only  waves 
propagating  outward  from  the  cylinder.  Therefore,  the  solution  to  the 
above  differential  equation  is  the  Hankel  Function  of  the  first  kind;  thus 

CO 

W-.  [2  2 

I’Q(r,<|);k)  •Y  -  k  )  r  cos  ncj)  (1^) 


The  constant  A  can  be  determined  from  the  boundary  condtion  at  the 
n 

fluid-shell  interface  where  continuity  requires  the  shell  and  fluid  par¬ 
ticle  velocity  to  be  equal,  l.e.. 


—  (((),z,t)  =  V  (r,(|),z,t)| 

*  r=a 


The  fluid  particle  velocity  v^  Is  related  to  the  fluid  pressure  by 


the  momentum  balance  relation 
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3Vf 
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Combining  Equations  (15)  and  (16)  while  suppressing  the  harmonic  time 
dependence  we  have 
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where 
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In  the  wave  number  domain  the  above  equation  can  be  expressed  as 
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where 


w  ((t);k)  =  y  cos  n4)  W  (k) 
o  n 
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Combining  Equations  (14)  and  (18)  and  solving  for  A  one  obtains 
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Consequently 
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where 
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is  the  transform  of  the  fluid  impedance.  The  transform  of  the  fluid  pres¬ 
sure,  p(a,(j),z)  acting  on  the  shell  surface  is  obtained  by  setting  r=a. 

Thus 
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DETERMINATION  OF  CYLINDER'S  RADIAL  RESPONSE 

Substituting  the  transformed  expressions  for  the  exciting  force  and 
fluid  pressure  acting  on  the  cylinder  into  Equation  (6c)  yields 
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The  radial  component  W  (k)  can  now  be  determined  from  Equations  (6a), 

n 

(6b)  and  (21)  by  applying  Cramer's  rule,  l.e.. 
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Z  +  Z  describes  the  Impedance  of  the  entire  cylindrical  shell, 
np  nm 

Note  that  Z^^  approaches  the  Impedance  of  a  plate  as  fi  approaches  “>  . 

On  the  other  hand,  Z  ,  the  Impedance  due  to  the  membrane  terms  of  the 

nm 

cylindrical  shell,  approaches  zero  when  a»l. 

By  applying  the  Inverse  Fourier  transform, 

+  OO 

W^(z)-^  Jw^(k)e^''^dk  (23) 

-  OD 

to  Equation  (19)  and  substituting  the  result  Into  Equation  (3c),  the 
radial  velocity  of  the  cylinder  can  be  expressed  as 
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w(<t>,Z,t)  =  - 2~ 

(27r)  a 


cos  n(()e  J  G^(k)dk 


(24a) 


where 


Ikz  ,, 

^  Z  (k)  +  Z  (k)  +  Z  (k) 
np  nm  na 


(24b) 


Similar  expressions  can  be  obtained  for  the  displacements  u  and  v. 


RADIATED  PRESSURE  FIELD 

The  radiated  pressure  field  can  now  be  expressed  as 


p(r,<t),2,t)  =  — -  ^  e  cos  n(()e  Z  (k)G  (k)dk 

n  nr  n 

^  n-0 


where  Z  (k)  and  G  (k)  are  defined  by  Equations  (19b)  and  (24b),  respec- 
nr  n 

tlvely . 

Referring  to  the  expressions  for  Z  (or  Z  )  it  can  be  seen  that 

nr  na 

for  k  >k,Z  has  both  a  real  (resistive)  and  imaginary  (reactive)  part, 
o  nr 

Thus  for  energy  will  be  radiated  into  the  acoustic  medium.  On  the 

other  hand,  when  k>k  the  arguments  of  Z  become  Imaginary  and  the  Hankel 

o  nr 

function  can  be  replaced  by  the  modified  Hankel  function  with  a  real 
argument,  namely. 


H  (lx)  -  -  K  (x) 

n  n 
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(27) 


Therefore  for  k>k 

o 


which  is  purely  a  reactive  and  mass-like  impedance.  Thus  for  k>k^  no 
energy  will  be  radiated. 

METHOD  OF  SOLUTION 

The  Integral  in  Equation  (24a)  may  be  referred  to  as  a  modal  Influence 
function  which  relates  the  response  of  the  cylinder  at  a  point  z  due  to  a 
force  applied  in  this  case  at  z=0  for  a  particular  mode  n.  Similarly  the 
integral  in  Equation  (25)  relates  the  radiated  pressure  to  the  force  for 
a  particular  mode  n.  The  total  response  or  radiated  pressure  is  then 
obtained  by  summing  the  modal  functions  for  all  the  circumferential  modes. 

There  is  no  known  closed  form  solution  to  the  above  integrals.  The 
integrands  are  transcendental  functions  further  complicated  by  the  presence 
of  poles  and  branch  points  in  the  complex  k  plane.  The  Integrand  is  not 
analytic  for  all  values  of  k  in  the  complex  k  plane.  The  technique 
employed  to  evaluate  the  modal  functions  is  similar  to  the  one  developed 
for  evaluating  the  Green’s  function  for  a  line  driven,  fluid  loaded  infinite 
plate.  This  method  involves  numerically  evaluating  the  Inverse  Fourier 
transform  Integrals  along  the  real  axis  after  structural  damping  has  been 
introduced  into  the  Integrand.  Structural  damping  is  introduced  by  replac¬ 
ing  Young's  modulus  E  by  E(l-lri)  and  is  necessary  to  move  the  singularities 
off  the  real  axis,  n  is  the  loss  factor  for  the  shell  material.  Briefly, 
the  numerical  scheme  is  based  on  using  variable  Integration  Increments  along 


the  path  of  Integration.  The  size  of  the  Increment  Is  automatically 
adjusted  so  as  to  keep  the  Incremental  change  of  the  integrand  values 
within  prescribed  limits.  In  essence,  the  second  derivative  of  the  function 
controls  the  mesh  size.  Thus  where  the  integrand  is  not  well  behaved  near 
a  singularity  the  mesh  size  is  very  fine.  On  the  band  where  the  integrand 
is  well  behaved  the  mesh  size  is  much  larger.  This  method  produces  a  set 
containing  a  minimum  number  of  integration  points  which  accurately  describes 
the  function  to  be  integrated  in  a  very  smooth  and  continuous  fashion. 

Once  the  set  of  points  has  been  generated,  the  integral  is  evaluated  by 
using  Simpson's  rule. 

DISCUSSION  OF  RESULTS 

The  results  presented  in  this  report  are  for  a  shell  thickness  to 
radius  ratio  of  h/a=.01.  Response  calculations  were  performed  for  both  the 
in  vacuo  and  water  loaded  cases.  In  order  to  put  the  results  into  perspec¬ 
tive,  they  are  compared  with  those  of  an  equivalently  loaded  plate.  This 
is  particularly  important,  since  an  analysis  of  this  type  for  a  point 
excited  cylindrical  shell  has  apparently  never  been  previously  published. 

MODAL  RESPONSE  OF  CYLINDRICAL  SHELL 

The  total  response  of  a  cylinder  is  a  summation  over  all  its  modal 
responses.  In  practice,  however,  only  a  finite  number  of  modes  have  to  be 
summed.  The  number  of  modes  to  be  summed  depends  primarily  on  the  excita¬ 
tion  frequency  and  the  degree  of  accuracy  desired.  Consequently,  it  is 
informative  to  see  how  individual  circumferential  modes  vary  as  a  function 
of  frequency  and  also  the  relative  contribution  of  the  various  modes  when 
the  cylinder  is  excited  at  a  particular  frequency. 


The  results  presented  In  this  section  all  have  been  normalized  to  an 
equivalent  line  loaded  plate.  The  drive  point  velocity  of  a  line  loaded 
Infinite  plate  In  vacuo  has  an  exact  solution  of  the  form 


G  (o) 
P 


(1+in)^^^ 


(1+1) 


^  (i+in/4)(i+i) 


(28) 


The  normalizing  factor  is  F'k  /Auan  where  F'  =  F/2iTa  has  the  units  of  force 

P 

per  unit  length  and  F  is  the  point  force  acting  on  the  cylinder. 

In  Figure  2,  the  real  and  imaginary  components  of  the  admittance  in 

vacuo  for  the  n=0  circumferential  mode  are  shown  as  a  function  of  the 

frequency  parameter  a.  The  frequency  parameter  ct  is  often  used  for  plates 

and  is  equal  to  k  /k  where  k  is  the  acoustic  wave  number  of  the  acoustic 
op  o 


where  is  the  coincident  frequency  of  the  plate.  In  terms  of  the  non- 

1/2 

dimenslonallzed  frequency  parameter  a  =  (Cp/c)(n8) 

Referring  to  Figure  2  the  maximum  response  as  a  function  of  frequency, 

a  resonance,  for  the  n=0  occurs  at  the  ring  frequency  of  the  cylindrical 

shell.  That  is  when  the  compresslonal  structureborne  wave  length  in  the 

shell  material  equals  the  circumference  of  the  cylinder.  In  vacuo  the  ring 

frequency  is  given  by  f  =  c  /2'n'a.  Above  the  ring  frequency,  the  admittance 

o  p 

converges  rapidly  to  the  admittance  of  a  line  excited  plate  which  has  a 
normalized  value  of  .975  +  11.025  with  a  loss  factor  of  .1.  For  n"0  the 
tangential  component  of  motion  for  the  shell  is  zero. 

In  Figure  3,  the  real  and  imaginary  components  of  the  admittance  of 


the  n-0  mode  and  a  line  excited  plate  are  compared  with  water  loading 
Included.  As  In  the  In  vacuo  case,  the  cylindrical  shell  admittance  differs 
considerably  from  that  of  a  line  driven  plate  in  the  vicinity  of  the  ring 


frequency  of  the  cylinder.  However,  the  amplitude  at  resonance  Is  much 
lower  with  fluid  loading  than  the  corresponding  In  vacuo  case.  This  Is  due 
to  the  fact  that  at  the  ring  frequency  the  structural  wave  length  Is  larger 
than  the  acoustic  wave  length  In  the  fluid  and  consequently  the  n=0  mode 
of  a  cylinder  Is  a  very  efficient  radiator  at  this  frequency  and  Is  strongly 
damped  by  radiation. 

In  addition,  the  presence  of  the  fluid  lowers  the  nondlmenslonalized 
frequency  of  the  n=0  mode  from  f2=l  for  the  In  vacuo  case  to  Q-.Tl  with 
fluid  loading.  Junger  and  Felt  have  obtained  an  expression  for  approxi¬ 
mating  the  natural  frequency  of  the  n=0  mode  of  a  submerged  cylindrical 
shell. ^  The  approximation  is 


It  is  reassuring  to  note  that  the  above  expressioi.  predicts  that  the 
reduced  ring  frequency  of  the  n=0  mode  in  water  to  be  77. 

In  Figure  4  the  admittance  corresponding  to  circumferential  mode 
number  6  is  shown  as  a  function  of  frequency  with  and  withov  t  water  loading. 
Unlike  the  n=0  circumferential  mode  there  are  now  two  resonant  frequencies. 
As  a  matter  of  fact  for  each  n>l  there  are  always  two  resonant  frequencies 
present  in  the  radial  response.  The  lower  resonant  frequency  Is  primarily 
an  Inextenslonal  type  mode  where  as  the  higher  resonant  frequency  is 
primarily  an  extenslonal  or  membrane  type  mode. 

In  Figure  5  the  drive  point  response  of  the  forty  first  circumferential 
mode  Is  Illustrated  as  a  function  of  frequency  with  water  loading  Included. 
Again  there  are  two  resonant  frequencies.  However,  unlike  the  case  for  low 


values  of  n,  there  is  very  little  radial  motion  due  to  extentlonal  or 
membrane  strains  at  the  higher  resonant  frequency.  This  difference  can 
be  explained  by  energy  considerations.^  Arnold  and  Warburton  have  shown 
that  at  low  circumferential  wave  numbers  n/a,  the  strain  energy  associated 
with  bending  is  low  compared  to  the  strain  energy  due  to  stretching  of  the 
reference  surface.  On  the  other  hand,  at  higher  circumferential  wave 
numbers  the  bending  strain  energy  is  high  whereas  the  stretching  strain 
energy  is  low. 

Since  the  response  of  the  cylinder  depends  upon  all  the  circumfer¬ 
ential  modes  of  vibration,  it  is  of  interest  to  see  how  each  individual 
mode  contributes  to  the  total  response  when  the  cylinder  is  excited  at  a 
particular  frequency.  Figure  6  demonstrates  the  relative  contribution  of 
the  various  modes  when  the  cylinder  is  driven  at  a  frequency  near  the 
resonance  frequency  corresponding  approximately  to  the  forty  first  bending 
mode.  As  can  be  expected,  the  modes  in  the  vicinity  of  the  resonant 
bending  mode  contribute  substantially  more  to  the  radial  response  than  do 
modes  farther  away.  The  radial  response  of  the  cylinder  at  this  exciting 
frequency  would  be  primarily  due  to  circumferential  bending  rather  than 
stretching  of  the  reference  surface.  It  can  also  be  shown  that  for  values 
of  n  somewhat  larger  than  the  predominantly  excited  modes,  the  admittance 
decreases  at  the  rate  of  1/n  (see  Appendix  B) .  Thus  it  is  relatively  easy 
to  minimize  truncation  errors  when  calculating  the  drive  point  admittance. 

RESPONSE  OF  A  POINT  EXCITED  CYLINDER 

As  previously  pointed  out  the  response  of  a  point  excited  cylinder 
should  converge  at  higher  frequencies  to  that  for  a  point  excited  plate. 
Thus,  for  convenience,  the  results  in  this  section  are  normalized  to  the 
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admittance  of  a  point  driven  plate  in  vacuo  (see  Appendix  A) .  In  vacuo, 


the  velocity  of  a  point  excited  plate  at  radius  r  can  be  expressed  as 

+  oo 


„  fU  (rk)dk 
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J  P 


(29) 


and  with  fluid  loading,  the  velocity  is 

+  oo 
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kJ  (rk)dk 
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Z  (k)  +  Z  (k) 
P  a 


(30) 


where  =  cylindrical  Bessel  function  of  zero  order 

Z  (k)  =  is  the  transform  of  the  plate  surface  impedance  operator 
P 

Z^(k)  =  is  the  transform  of  the  fluid  Impedance 
r  =  distance  from  the  applied  force 
It  can  be  shown  that  the  admittance  or  drive  point  velocity  (r=0) 
of  the  plate  in  vacuo  is  Independent  of  frequency  and  can  be  expressed  as 


V  (0) 
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^  p  c  h^ 
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(31) 


The  above  term  is  used  to  normalize  all  results  presented  in  this 
section.  The  velocity  response  away  from  the  drive  point  for  the  in  vacuo 
and  fluid  loading  cases  are  obtained  by  numerically  integrating  Equations 
(29)  and  (30),  respectively,  by  a  method  similar  to  the  one  outlined  for  a 
cylinder.  In  the  case  of  fluid  loading  the  admittance  must  also  be 
obtained  numerically. 

In  Figure  7  the  normalized  admittance  for  a  point  excited  cylinder 
is  shown  with  and  without  water  loading  as  a  function  of  the  frequency 
parameter  ot.  These  results  are  compared  with  the  admittance  of  a  point 


driven  plate.  Since  the  admittance  of  a  plate  in  vacuo  is  independent  of 
frequency,  its  normalized  value  is  therefore  a  constant  equal  to  one.  As 
should  be  expected,  the  greatest  difference  between  a  plate  and  cylindrical 
shell  occurs  in  the  vicinity  of  the  ring  frequency  and  at  lower  frequencies. 
However,  it  is  surprising  how  quickly  the  admittance  of  the  cylinder  con¬ 
verges  to  a  plate  for  both  the  in  vacuo  and  fluid  loaded  cases.  Neverthe¬ 
less,  the  frequency  at  which  the  admittance  of  a  point  driven  cylinder 
approaches  that  of  a  plate  depends  on  the  parameter  h/a.  For  values  of 
h/a  larger  than  used  here,  the  frequency  would  be  higher.  On  the  other 
hand,  for  smaller  values  of  h/a  the  convergence  frequency  would  be  lower. 

The  peaks  in  the  admittance  curves  occur  at  frequencies  corresponding 
to  the  ring  frequency  of  the  cylinder.  The  response  at  this  frequency  is 
dominated  by  the  n=0  mode.  At  lower  frequencies,  i.e.,  low  values  of  a, 
the  admittance  curves  for  the  cylinder  are  shown  as  smooth  functions  of  a 
for  convenience  only.  To  further  explain  this,  let  us  assume  that  the 
response  of  the  cylinder  has  no  z  dependence;  then  the  natural  frequencies 
of  the  circumferential  bending  modes  in  vacuo  are  given  by 


f 
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^P  Bn  (n^+1) 

2n  a  ,2  ,.1/2 
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(32) 


The  above  relation  can  also  be  used  to  approximate  exciting  frequen¬ 
cies  which  are  half  way  between  any  two  adjacent  circumferential  bending 
modes  of  the  cylinder.  Figure  8  Illustrates  the  effect  on  the  velocity 
response  when  the  exciting  frequency  coincides  with  the  frequency  of  a 
resonant  bending  mode  (n=lnteger)  and  when  it  does  not  (for  instance  when 
n=lnteger  +  .5).  As  can  be  seen  from  the  figure  there  is  an  oscillation 
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in  the  response  at  lower  frequencies.  However,  as  the  frequency  increases, 
the  disparity  becomes  less  and  less  until  the  response  curve  becomes  smooth 
regardless  of  the  exciting  frequency.  This  is  because  the  circumferential 
bending  waves  are  dispersive  and  their  wave  lengths  become  smaller  as  the 
exciting  frequency  increases.  Consequently,  at  higher  frequencies  the 
waves  damp  out  more  quickly,  thus  eliminating  any  reinforcement  or  cancel¬ 
lation  which  would  tend  to  occur  at  lower  frequencies.  It  should  be 
pointed  out,  however,  that  the  rate  at  which  the  envelope  closes  depends 
on  the  damping  value. 

Figures  9  and  10  demonstrate  the  differences  in  the  velocity  profile 
away  from  the  drive  point  between  a  plate  and  cylinder  when  excited  at  the 
ring  frequency  of  a  cylinder  without  and  with  water  loading,  respectively. 
For  the  cylinder  the  velocity  profile  is  along  the  z  axis  and  the  distance 
from  the  drive  point  is  expressed  in  nondimens Iona 11 zed  form  as  z/a. 

Figures  9  and  10  show  considerable  differences  in  the  velocity  profiles 
between  a  point  driven  plate  and  cylinder,  particularly  for  the  in  vacuo 
case.  This  should  be  expected  since,  at  the  ring  frequency,  the  radial 
motion  of  the  cylinder  is  primarily  due  to  extensional  strains  of  the 
shell's  middle  surface.  However,  water  loading  does  moderate  the  differ¬ 
ences  in  the  velocity  profile  between  a  point  driven  plate  and  a  cylinder. 
In  any  event,  as  the  admittance  of  a  point  driven  plate  and  cylinder  ap¬ 
proach  each  other  at  higher  frequencies,  one  would  also  expect  their 
respective  velocity  profiles  in  the  vicinity  of  the  drive  point  to  become 
increasingly  similar.  The  degree  of  similarity  will  depend  not  only  on  the 
exciting  frequency  but  also  on  the  amount  of  damping,  for  it  is  the  damping 
value  which  determines  how  quickly  the  propagating  waves  from  the  drive 
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point  decay  and  therefore  how  localized  the  vibration  pattern  becomes.  The 
more  localized  the  vibration  pattern  is,  the  less  is  the  effect  the  shell 
curvature  has  upon  the  vibration  response. 

Figure  11  Illustrates  the  resulting  velocity  profile  in  vacuo  around 
the  circumference  of  the  cylinder  when  the  excitation  frequency  coincides 
with  the  ring  frequency  of  the  cylinder.  At  the  ring  frequency,  the 
response  of  the  cylinder  is  a  combination  of  the  predominantly  excited 
bending  modes  superimposed  upon  the  predominantly  excited  membrane  modes. 
Figure  12  illustrates  the  contribution  of  the  various  mtdes  to  the  total 
drive  point  response  at  the  ring  frequency.  The  resonant  peak  that  occurs 
at  the  lower  mode  numbers  is  due  to  the  resonant  membrane  modes  as  well  as 
the  rigid  body  mode  (n=l) .  On  the  other  hand,  the  second  peak  is  due  to  a 
resonant  bending  mode.  The  dominant  resonant  bending  mode  has  an  n  value 
of  18.  This  mode  has  a  structural  wavelength  that  corresponds  to  an  angle 
of  20  degrees  along  the  circumference  of  the  cylinder.  As  a  result,  the 
lobes  in  the  response  occur  at  approximately  20°  intervals.  If  fluid 
loading  is  included  (Figure  13),  the  response  at  the  ring  frequency  is  not 
only  substantially  reduced  but  the  lobes  rapidly  disappear  away  from  the 
drive  point.  This  can  be  attributed  to  the  fact  that  in  fluid  the  excited 
membrane  modes  of  the  cylinder  are  critically  damped  due  to  radiation.  To 
further  illustrate  this,  the  contribution  of  each  mode  to  the  total  drive 
point  response  with  fluid  loading  is  shown  in  Figure  lA.  A  comparison 
with  the  in  vacuo  case.  Figure  12,  clearly  demonstrates  that  the  presence 
of  the  fluid  greatly  reduces  the  amplitude  of  the  resonant  peak  due  to  the 
resonant  membrane  modes  and  the  rigid  body  mode.  This  reduction  takes 
place  because  these  modes  are  efficient  radiators  and  therefore  lose  energy 
through  radiation  damping. 
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Figures  15  and  16  compare  the  velocity  response  of  a  cylinder  and 


plate  away  from  the  drive  point  along  the  generator  without  and  with 
water  loading,  respectively,  at  an  exciting  frequency  corresponding  to 
a=.5.  As  the  two  figures  indicate,  at  a=.5  the  response  of  a  plate  and 
cylinder  are  identical  at  least  for  the  range  of  z/a  where  calculations 
were  performed.  Thus,  for  the  values  of  h/a  and  damping  used  here,  the 
vibration  response  of  a  cylinder  can,  for  all  practical  purposes,  be  ac¬ 
curately  represented  by  a  flat  plate  in  the  vicinity  of  a=.5. 

Figures  15  and  16  also  illustrate  the  velocity  response  of  a  point 
driven  plate  excited  at  a  frequency  of  a=1.5  with  and  without  water 
loading,  respectively.  In  view  of  the  results  for  a=.5,  the  velocity 
response  of  a  point  driven  cylinder  should  be  and  in  fact  is  identical  to 
a  point  driven  plate  for  a=1.5. 

So  far  only  the  velocity  response  along  the  z  axis  of  the  cylinder  has 
been  compared  with  that  of  a  point  driven  plate  of  equivalent  radius  r. 
Indeed,  if  at  some  frequency  a  cylinder  can  be  approximated  by  a  plate, 
the  velocity  response  along  the  circumference  of  a  cylinder  away  from  the 
drive  point  should  also  be  similar  to  a  point  driven  plate.  In  Figure  17 
such  a  comparison  is  made  between  a  point  driven  cylinder  and  plate  with 
a=.5.  From  Figure  17  it  is  obvious  that  the  two  responses  are  indeed  very 
similar  for  the  range  of  z/a  considered  and  for  a  damping  value  of  .1. 

Figure  18  demonstrates  how  localized  the  velocity  pattern  near  the 
drive  point  becomes  with  a=.5.  In  Figure  18  a  10°  increment  corresponds 
to  a  z/a  increment  of  approximately  .175  in  Figure  17. 

All  the  previous  results  have  been  arrived  at  by  assuming  a  damping 
factor  .1.  Although  a  damping  factor  of  .1  may  be  a  reasonable  value 
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for  many  structures.  It  may  be  too  high  for  others.  Consequently,  in  order 
to  demonstrate  the  effect  of  structural  damping  on  the  vibration  response 
of  a  cylinder,  calculations  were  also  performed  using  a  very  small  damping 
value  of  .01.  In  Figure  18  the  velocity  envelope  around  the  circumference 
of  a  cylinder  is  shown  when  the  excitation  frequency  corresponds  to  the 
ring  frequency  In  vacuo.  A  comparison  between  Figures  11  and  19  demon¬ 
strates  how  the  velocity  envelope  becomes  more  complicated  when  the  damping 
value  Is  changed  from  .1  to  .01.  This  Is  due  to  the  fact  that  structural 
waves  generated  at  the  drive  point  do  not  damp  out  as  rapidly  and  thus 
travel  around  the  circumference  of  the  cylinder  In  opposite  directions  with 
larger  amplitudes.  Consequently  when  the  two  waves  meet,  they  either  rein¬ 
force  or  interfere  with  each  other  depending  on  the  observation  point  on 
the  circumference  of  the  cylinder.  As  previously  pointed  out,  the  pre¬ 
dominantly  excited  bending  mode  at  the  ring  frequency  in  vacuo  Is  n=18  and 
therefore  has  a  structural  wavelength  corresponding  to  20  degrees  along  the 
circumferences  of  the  cylinder.  Therefore,  in  the  vicinity  of  (J)“180°  where 
the  two  traveling  waves  have  nearly  the  same  magnitude,  the  velocity 
envelope  has  a  periodic  wave  pattern  but  with  a  wavelength  of  10°  along 
the  circumference  due  to  the  Interaction  of  the  two  waves  traveling  in 
opposite  directions  around  the  cylinder.  Closer  to  the  drive  point  the 
velocity  profile  becomes  more  complicated  because  of  the  greater  contri¬ 
bution  of  the  membrane  modes  to  the  total  response.  Figure  19  also  shows 
the  actual  traveling  wave  pattern  at  a  time  "T"  when  the  drive  point 
response  Is  at  Its  maximum  amplitude. 

At  frequencies  above  the  ring  frequency  where  the  membrane  modes  have 
little  or  no  effect  upon  the  response,  the  velocity  profile  is  less 


complicated  for  a  damping  value  of  .01.  This  fact  Is  Illustrated  In 
Figure  20  where  the  cylindrical  shell  In  vacuo  Is  point  driven  at  a  fre¬ 
quency  corresponding  to  a=.5.  For  convenience,  the  maxima  and  minimum 
points  of  the  velocity  response  are  represented  by  an  envelope  as  a  function 
of  the  circumferential  coordinate  (j>. 

As  the  observation  point  moves  further  away  from  the  drive  point  the 
envelope  widens.  The  maximum  width  of  the  envelope  would  occur  at  (J)=180 
where  the  two  traveling  waves  are  of  equal  magnitude.  However,  closer  to 
the  drive  point  the  standing  wave  pattern  of  the  velocity  envelope  Is  less 
apparent  since  the  response  Is  primarily  caused  by  the  wave  traveling  the 
least  distance  from  the  drive  point.  For  comparison  purposes  the  response 
of  a  point  excited  plate  with  a  structural  damping  value  of  .01  would  tend 
to  predict  an  average  value  for  the  cylinder's  response  In  Figure  20. 

Away  from  the  drive  point  along  the  generator  of  the  cylinder,  Inter¬ 
ference  of  the  waves  traversing  the  cylinder  also  occurs.  However,  the 
effect  of  this  interference  on  the  response  of  the  cylinder  Is  always  less 
than  for  a  case  where  the  waves  circumscribe  the  cylinder  at  the  drive  point. 
This  is  because  along  the  generator  away  from  the  drive  point,  the  waves 
must  follow  a  helical  path  to  the  observation  point.  The  distance  of  the 
length  of  the  helical  path  is  always  larger  than  2TTa,  the  circumference  of 
the  cylinder.  Consequently,  due  to  damping,  the  effect  of  interference  at 
the  observation  point  Is  diminished. 

SUMMARY 

In  this  report  an  analysis  has  been  presented  for  calculating  the 
vibration  response  of  a  point  excited.  Infinitely  long,  thin  cylindrical 
shell  Immersed  in  an  acoustic  medium.  Using  the  Fourier  transform  approach 


and  refined  numerical  techniques  it  was  possible  to  obtain  a  direct  solu¬ 
tion  without  the  Inhibiting  assumptions  or  approximations  used  In  may  pre¬ 
vious  studies.  In  order  to  verify  the  analysis  to  the  best  possible  extent, 
results  were  compared  with  those  of  a  fluid  loaded,  point  excited  plate  at 
higher  frequencies.  The  vibration  response  of  the  point  driven  plate  was 
obtained  by  employing  a  numerical  technique  similar  to  the  one  used  to 
obtain  a  solution  for  the  cylinder. 

Example  calculations  were  performed  for  a  cylinder  with  a  wall  thick¬ 
ness  to  radius  ratio  of  .01.  A  comparison  between  the  results  for  a  point 
excited  cylinder  and  plate  has  given  additional  Insight  Into  the  frequency 
regions  where  a  cylinder  can  be  reasonably  approximated  by  a  plate  (and 
where  it  cannot  be).  It  was  found,  not  unexpectedly,  that  the  greatest 
difference  between  the  responses  of  a  point  excited  plate  and  cylinder 
occurs  at  low  frequencies  and  also  in  the  vicinity  of  the  ring  frequency  of 
the  cylinder  where  membrane  type  modes  dominate.  However,  It  was  surprising 
how  quickly  the  two  responses  approach  each  other  above  the  ring  frequency 
but  well  below  the  critical  frequency  of  the  plate  (a=1.0).  Consequently 
for  the  example  presented,  that  Is  with  an  h/a  ratio  of  .01,  the  response 
of  a  point  excited  cylinder  can  be  approximated  with  reasonable  accuracy  by 
the  response  of  a  point  excited  plate  in  the  mid  frequency  range  (a=.5) 
as  well  as  at  higher  frequencies  (of  evaluation) .  It  is  very  advantageous 
whenever  the  response  of  a  cylinder  can  be  approximated  by  a  plate  since 
the  analysis  of  a  plate  Is  much  easier  and  cheaper  than  for  a  cylinder. 

Although  the  radiated  pressure  in  the  near  field  for  the  cylinder  was 
not  calculated.  It  can  be  done  using  Equation  (25).  The  numerical  method 
used  to  evaluate  the  pressure  Integrals  would  be  Identical  to  the  one 


developed  to  evaluate  the  response  integrals.  It  would  be  interesting  to 
compare  the  radiated  pressure  field  of  a  point  excited  cylinder  and  plate, 
particularly  in  the  frequency  region  where  the  membrane  or  extentlonal 
strains  in  the  cylinder  are  of  significance. 


APPENDIX  A 


NORMALIZED  RESPONSE  OF  A  POINT  EXCITED  CYLINDRICAL  SHELL 

The  results  presented  in  this  report  have  been  normalized  in  two 
different  ways.  First  of  all,  when  investigating  the  behavior  of 
individual  circumferential  modes,  the  response  was  normalized  to  the 
admittance  of  a  line  excited  plate  in  vacuo.  On  the  other  hand,  the 
total  response,  which  contains  all  the  circumferential  modes  of  vibra~ 
tion,  was  normalized  to  the  admittance  of  a  point  excited  plate  in  vacuo. 


NORMALIZING  TO  A  POINT  EXCITED  PLATE 

The  expression,  Equation  (21a),  for  the  response  of  a  point  excited 
cylindrical  shell  can  be  written  in  the  follow  form 
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The  drive  point  velocity  of  a  point  excited  plate  in  vacuo  is 


v  (0) 


VT_F 

4 


p  c  h 
s  p 


(A2) 


27 


Dividing  Equation  (Al)  by  Equation  (A2)  the  normalized  cylinder  response 
may  be  expressed  as 
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NORMALIZING  TO  A  LINE  EXCITED  PLATE 

The  drive  point  velocity  of  a  line  excited  plate  in  vacuo  without 
damping  has  an  exact  solution  of  the  form 
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If  Equation  (Al)  is  divided  through  by  2^— ^ ,  the  normalized  admit 
tance  of  the  nth  mode  for  a  point  excited  cylinder  is  given  by 

+  00 


w  (0,0, t)  =  |(il)  (f26)^^^  G 
n  "  \a  /  n 


r  dk 
I  Z  (k)  + 

J  ns 


Z  (k) 
ns  na 


(A5) 


where  F’ 


ZTTa 


Similarily,  the  normalized  admittance  for  a  line  excited  plate  in 


vacuo  IS 


APPENDIX  B 


HIGH  FREQUENCY  APPROXIMATION  FOR  A  POINT  EXCITED  INFINITELY  LONG  CYLINDER 
IN  VACUO 

Heckl  has  obtained  a  high  frequency  approximation  for  the  admittance 
of  a  point  excited  thin  cylindrical  shell  In  vacuo. ^  For  a  particular 
circumferential  mode  Heckl 's  expression  Is  given  by 
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The  above  equation  Is  presumably  valid  in  the  frequency  region  where 
Q»l.  The  admittance  or  drive  point  velocity  of  the  shell  can  be 
expressed  as 
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n=0 

In  Figure  21  the  high  frequency  approximation  is  compared  with  the 
admittance  previously  obtained  for  a  point  excited  cylinder  and  plate. 

It  is  evident  that  Heckl 's  high  frequency  approximation  converges  to  a 
point  excited  plate  and  is  an  excellent  approximation  to  a  cylinder  in  the 
frequency  range  of  a>.5  for  a  h/a  ratio  of  .01.  It  is  also  interesting 
to  note  that  for  oi=.5,  ^=.7  which  indicates  the  above  approximation  is 
valid  for  frequencies  much  lower  than  originally  assumed.  However,  it 
is  apparent  that  Heckl 's  approximation  neglects  the  membrane  modes  of  the 
cylinder. 
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NORMALIZED  VELOCITY 


Figure  2  -  Admittance  of  the  Circumferential  n=0  Mode 
for  a  Cylindrical  Shell  in  Vacuo 
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Figure  4  -  Admittance  of  a  Cylindrical  Shell  Immersed  in  Water  and 
in  a  Vacuum  for  Circumferential  Mode  6 
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Figure  7  -  Drive  Point  Velocity  of  a  Cylindrical  Shell  and  a  Plate  as  a  Function  of  Frequency 


Comparison  of  the  Velocity  Responses  for  a  Cylinder  and  Plate  In  Vacuo,  Point 
Excited  at  the  Ring  Frequency 


Comparison  of  the  Velocity  Responses  for  a  Cylinder  and  Plate  in  Vacuo,  Point 
Excited  at  Frequencies  of  a=0.5  and  a=1.5 
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Figure  19  -  Normalized  Velocity  Profile  Along  the  Circumference  of  a 
Cylindrical  Shell  with  a  Structural  Damping  Value  of  0.01 
and  Point  Driven  at  the  Ring  Frequency  in  Vacuo 


Figure  20  -  Envelope  of  the  Minima  and  Maxima  Response  of  a  Cylindrical  Shell  in  Vacuo  as 

Function  of  the  Circumferential  Coordinate  (}) 
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Figure  21  -  Drive  Point  Velocity  of  a  Cylindrical  Shell  and  a  Plate  in  Vacuo 
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